Ceo in intense femtosecond laser pulses: nonlinear dipole response and ionization 
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We study the interaction of strong femtosecond laser pulses with the Ceo molecule employing time- 
' dependent density functional theory with the ionic background treated in a jellium approximation. 

The laser intensities considered are below the threshold of strong fragmentation but too high for 
perturbative treatments such as linear response. The nonlinear response of the model to excitations 
by short pulses of frequencies up to 45 eV is presented and analyzed with the help of Kohn-Sham 
orbital resolved dipole spectra. In femtosecond laser pulses of 800 nm wavelength ionization is found 
to occur multiphoton-like rather than via excitation of a "giant" resonance. 

PACS numbers: 36.40.-c, 33.80. Rv, 31.15.Ew 
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I. INTRODUCTION 



Intense laser atom interaction exhibits nonlinear phenomena such as above threshold ionization (ATI), high harmonic 
• generation (HHG), and nonsequential multiple ionization (NSI) (see ^ for recent reviews). While some of those 
. ^ ' features are accessible in terms of a sequential "single active electron" (SAE) approach others are clear manifestations 
^ of many electron effects and correlation, e.g., NSI. The full ah initio solution of the time-dependent Schrodinger 
equation (TDSE) for two active electrons interacting in their full dimensionality with the laser and their parent ion 
^ ^ is already at the limit of what is possible with modern computers Q. Treating many electron systems in laser fields 
I— I thus needs further approximations. Density functional theory (DFT), extremely successful in electronic structure 
calculations of many-electron systems ( see, e.g., ^), has been extended to the time-dependent case (TDDFT) (see 
[Q, and, e.g., 1^ for a review). Despite the fact that TDDFT still lacks an equally solid foundation compared to that 
^ on which ground state DFT is built, it was successfully applied to metal clusters in laser pulses (see |^ for a review). 
Problems mainly arise when observables have to be extracted which are not obvious functionals of the Kohn-Sham 
\ orbitals or the total electron density, like in the study of NSI of atoms within TDDFT |^,^ , or when the results are very 
\^ • sensitive to the choice of the approximation to the unknown exchange-correlation potential. Compared to laser atom 
' interaction, in big molecules or clusters additional degrees of freedom are introduced: electronic degrees of freedom, 
including collective effects such as the formation of plasmons, vibrational degrees of freedom, or fragmentation. With 
laser pulses of different duration the equilibration of energy among the various channels can be probed. For Cgo this 
was nicely demonstrated in Ref. where the photoelectron spectra in fs laser pulses exhibited ATI peaks, a signature 
for direct multiphoton processes, which disappeared for longer pulses where collective effects set in. Concerning the 
^ \ ionization mechanism of Cgo in fs laser pulses there is a discrepancy in the literature. While in the recent work of 



Tchaplyguine et al. |10| from ion yield-curves vs. laser intensity direct multiphoton ionization was found to be the 
responsible pathway for ionization of Cgo, in an earlier publication Hunsche et al. pd| ] claimed it is the excitation of 
a giant resonance near 20 eV. Such a resonance at 20 eV in Cgo was first predicted theoretically by Bertsch et al. |T^ ] 
. ^ ^^'-^ confirmed later in an experiment by Hertel et al. p3| using synchrotron radiation for single-photon ionization 
K> , measurements. When compared to metal clusters where collective resonances occur at a few eV a 20 eV giant resonance 
5_] of about 12 eV width is quite remarkable. 

^ The nonlinear TDDFT treatment of Cgo in a laser pulse is numerically rather demanding because one has to allow 
for ionization which implies the use of a big numerical grid in order to represent the time-dependent Kohn-Sham (KS) 
orbitals properly. It is thus impossible, at least with computers nowadays available, to achieve both a detailed account 
of the soccer ball-like structure of Cgo and an accurate propagation of freed electron density far away from the ion 
and, possibly, back. Therefore we restrict ourselves to a jellium approach for the ionic background of the Ceo- Such a 
jellium model was employed in to study the photo absorption of atoms inside Cgo within linear response theory. 
It was found to share many of the relevant features with more demanding "first principle" calculations (like, e.g., in 
[ p2t ) and experiment . Jellium models were also successfully applied to metal clusters (see ||l^ for a review) . 
The paper is organized as follows. In Section O we present our TDDFT jellium model of laser Cgo interaction. In 



Section HI we characterize the dipole spectrum of our model in terms of single particle transitions. In Section IV 
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we study the dipole response and ionization after excitation by short pulses of different frequencies. In Section ^ we 
examine the ionization mechanism of our Cgo model in 800 nm fs laser pulses. Section ^ contains a brief summary 
and conclusion. 

II. THE MODEL: STATIC PROPERTIES 

The time-dependent Kohn-Sham (TDKS) equation for the orbital ^'icr(»",i) reads (atomic units are used unless 
noted otherwise) 

i^v|/^,(r,i)- (^-^V' + V{r)+yj{t) + V,,,[n^{r,t),n^{r,t)]^ *4.(r,t). (1) 

Here, a =1, j indicates the spin polarization, V{r) is the potential of the ions, V/(t) is the laser in dipole approximation, 
and Vcca ['^T > "il effective electron-electron interaction potential which is a functional of the electron spin densities 

n,(r,i) -^|*,,(r,t)|2. (2) 

4=1 

is the number of orbitals occupied by KS particles with spin c. The total electron density is 

nir,t) ^^n„{r,t). (3) 

The electron-electron potential is splitted, 

Vee^[ni,ni] = U[n] + V^cA^-i ,ni], (4) 



UM-ldV^ (5) 



where U[n] is the Hartree potential 



and V^ca ['^T ' "i] exchange correlation (xc)-part. Although the Runge-Gross theorem ensures that, in principle, 
the time-dependent KS scheme could yield the exact density n{r,t) on which all observables depend, in practice an 
approximation to the exchange-correlation potential V^ccr["-Ti '^il to be made. We chose the Slater expression 

t) = E ^^4^-xc4.(r, t), (6) 

where Wxcio- — ['^T'"'i] ~ C^l^itr] ~ [nicr,0], i.e., the self-interaction is removed, and the exchange-only 

local density approximation (XLDA) was employed. In TDDFT simulations of nonsequential ionization in laser atom 
interaction |^ we found that, compared to plain XLDA, the Slater potential (^) improved the ionization potentials and 
the fulfillment of Koopman's theorem significantly. Instead, for the kind of studies presented in this paper it is actually 
not necessary to go beyond XLDA. Inaccuracies in the, e.g., ionization potential for the outermost electron can be 
easily compensated by adjusting the free parameters of the jellium potential for the ionic background (see below). In 
Ref. Il^, for sodium clusters in laser pulses a self interaction correction was found to cause "minor differences" in 
dipole spectra leaving the ionization mechanism "essentially unchanged." Adding a correlation potential (we used the 
one proposed by Perdew and Wang |l^) had negligible effects on the results presented in this work. 

Usually the dipole approximation is said to be justified when the system under study is small with respect to 
the laser wavelength. However, when a light beam impinges on an object smaller than the wavelength the interior 
could be, nevertheless, field-free, provided the electron density is sufficiently high. This kind of screening would not 
be included in a TDDFT treatment with dipole approximation. The condition for laser light of wavelength and 
frequency cui being able to penetrate a plasma layer of thickness i5 is ^ = ■n(u}p/iOiY5/\i < 1 jlj]. This is the case for 
our parameters. Therefore, the dipole approximation is safe also in this respect. 

In our numerical code, the KS orbitals are expanded in spherical harmonics Yf^{{},(f). If the ground state has 
a closed shell structure the effective KS potential is spherical. Hence £ and m are "good" quantum numbers for 
the ground state configuration. In a linearly polarized laser field (in dipole approximation) the quantum numbers 
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m remain good, i.e., there is no m-coupling through the laser because the azimuthal symmetry is retained. V/(i) 
introduces an £± 1-couphng only. However, orbitals with different \m\ behave differently in the laser field. The radial 
KS wave functions are discretized in position space. Each time step, the effective potential has to be calculated which 
makes a TDKS solver significantly more time-consuming than an ordinary TDSE code running a corresponding SAE 
problem. The effective potential Vcc was expanded up to the dipole. Consequently, both the laser and V^c lead to 
an ^ ± f-couplingonly. We are confident that neglecting higher order multipoles of Vcc does not affect the validity of 
our conclusions p9[ | . The actual propa gat ion is performed in velocity gauge using an algorithm similar to the one for 
the TDSE described by MuUer in Ref . po|| . Probability density reaching the radial boundary of the numerical grid at 
~ 100 a.u. was removed by an imaginary potential. The eventually decreasing norm then can be interpreted as one 
minus the ionization probability of that orbital. 

The laser is polarized in z-direction so that in velocity gauge we have 

V,(t)^-iA(t)^ (7) 

where A{t) is the vector potential and the A'^-term has been transformed away (see, e.g., This potential leads 

to the above mentioned £ ±1 coupling. 

The ionic background is treated in a jellium approximation, i.e., the ions arc thought of being smeared out over 
a spherical shell with outer and inner radius Ro and Ri, respectively. The ionic charge density is constant for 
Ro > r > Ri and zero otherwise. The radii Ro, Ri are centered around the known radius of the Cgo-fuUcrcne, 
(i?o + Ri)/'^ = R = 6.7 a.u. In real Cgo there are 60 7r-electrons and 180 cr-electrons. Therefore the charge of the 
jellium background should be 240 a.u. However, 240 KS particles do not yield a self-consistent closed shell-structure 
for the ground state of our model. Since partially filled shells would spoil the spherical symmetry of the ground state 
we take 250 electrons instead (see also |l^]) which leads to a self-consistent closed shell ground state of the jellium 
model. Thus, introducing 

^«"' = ]^3^' ^ = 250, (8) 

and allowing for an additional potential depth parameter vq we arrive at an ionic potential 

r ^r-^3iRl~Rf)/2 for r < R,, 

Vir) = I -r;^ (3i?2/2 - [r72 + i?f /r] ) - i;o for R^ < r < Ro, 
{ -~r-^{Rl- Rl)/r for r > i?^. 

The parameters Ro, Ri, and vq can be varied in order to obtain a jellium-Cgo ground state which shares the relevant 
features with first-principle calculations of "real" Ceo (see also Ref. [Q). E.g., by decreasing Ro — Ri the overlap of 
the a and the tt single particle levels decreases. With vq the absolute position of the level scheme can be adjusted in 
order to meet, e.g., the ionization potential for the outermost electron. We chose Ro = 8.1, Ri = 5.3, and vq — 0.78. 
Some of the ground state properties of the model are shown in Fig. |l|. Because of the centrifugal barrier £{i + l)/2r^ 
the total potential is ^-dependent, and states with high £ are pushed outwards. The energy levels are 2{2£ + 1)- 
degenerated. The 250 KS particles can be subdivided in 200 n = 1-states (the cr-electrons) occupying ^-values from 
up to 9, and 50 n = 2-states (the 7r-electrons) . There are also bound n = 3-states ((5-levels) but they are not occupied 
in the ground state configuration. The orbital densities are also shown in Fig. |l]. Each 7r-electron wavefunction has a 
node near the jellium-shell radius R. The values of the single KS particle orbital eigenenergies are given in Table |. 
The highest occupied state is the tt state with £ = 4. From Koopman's theorem we therefore expect an ionization 
energy of — 0.274. Calculating the ionization energy by subtracting the total energy of C^q from that for neutral 
Ceo we obtain (on our numerical grid with a grid spacing Ar — 0.3, and the parameters Ri, Ro, vq as chosen above) 
/+ = 0.279 (7.59 eV) which agrees reasonably well with the former value and experiment (« 7.6 eV, [Q). In any case 
we expect for 800 nm laser light five photons being necessary for removing the outer electron. However, collective 
effects might occur so that more photons are required. In fact, there is an unresolved discrepancy in the literature 
about whether ionization of Ceo in the fs- regime works multiphoton-like |0 or through the excitation of a 20 eV giant 
resonance | [ll| . The results from our model concerning this question will be presented in Section 

In Table I we also enumerated the KS orbitals for the sake of easy reference later on. Since in each ^-shell there 
are £ + 1 different |r7i|-values we need 70 KS orbitals to describe our jellium-model interacting with the laser field. In 
each €-shell the KS orbitals are labeled from m = up to \m\ — £. Thus, e.g., orbital no. refers to the two electrons 
of opposite spin which, in the ground state configuration, populate the n = l, £ = m — cr-state whereas orbital no. 
69 is initially a pure n = 2, £ = \m\ = 4 7r-state, populated by four electrons with different spin and/or sign of m. 



3 



III. DIPOLE RESPONSE 



Within linear response theory excitations and giant resonances are commonly inferred from the photo absorption 
cross section, i.e., the imaginary part of the polarizability. Here we follow the different route of nonlinear TDDFT 
which allows us to distinguish between ionization, single particle transitions, or plasmons, and also accounts for higher 
order processes beyond single particle hole excitations. 

First we discuss the nonlinear response of our Cgo model to a kick at time i = by a delta-like electric field 
E{t) = AS{t). Such a kick is equivalent to the method proposed by Yabana and Bertsch in Ref. where the 
ground state wavefunctions are perturbed by giving them a coherent velocity field, i.e., — ^ exp(ifcz)^'j, in order to 
initiate a dipole moment evolving in time afterwards. From the Fourier-transform of the dipole d(t) = J zn{r) d^r 
the frequency dependent response is obtained. In Fig. ^ the Kohn-Sham orbital resolved dipole (KSORD) response for 
a delta-kick of magnitude A — 0.01 a.u. is presented. Such a kick can still be considered a small perturbation to the 
ground state configuration because only a fraction 3.2 x lO"'^ of the total electron density was freed (corresponding to 
an ionization probability of 8% for the outermost electron). On the left a contour plot of the dipole vs. orbital number 
and time is shown. On the right-hand-side the corresponding spectrum is plotted, obtained by Fourier-transforming 
the dipoles of the individual KS orbitals. The total dipole and its Fourier transform are also included in Fig.||. 
Looking at the dipoles vs. time one easily identifies the different ^-shells. The cr-electrons are KS orbital numbers 
0-54, the TT-electrons range from 55-69. Examining the KSORD spectra on the right one clearly identifies rather 
narrow vertical lines. Each of those vertical lines can be understood as a single particle transition between ground 
state KS levels. For a certain KS orbital the dipole strength is particularly high for those transitions where one of the 
two levels involved is the one which is occupied in the ground state configuration of this KS orbital. This explains 
the parabola-like structures of strong dipole emission visible in the KSORD spectra. 

The low frequency lines in the range 0-5 eV stem from transitions of the type I l±l with the ri-quantum number 
fixed. Lines between 5 and 15 eV are caused by transitions of type -kI —> S{£ ± 1), and are therefore particularly 
pronounced for orbitals representing 7r-electrons initially. The lines with high dipole strength along a parabola-like 
structure around 20 eV for both, tt- and cr-electrons originate from transitions of the type a£ 7t{£+1) oi iri a{£~l) 
(right branches) and ai — > Tr{£ — 1) or tt£ — + (t{£ + 1) (left branches). Even higher dipole emission around 30 eV stems 
from CT(5-transitions. A common feature for all orbitals where £ = \m\ initially (i.e., orbitals number 0, 2, 5, 9, 14, 20, 
27, 35, 44, 54, 55, 57, 60, 64, 69) is the relatively strong emission around 30 eV. The reason why the £ = |m|-orbitals 
preferably radiate at those frequencies is that £ ^ £ — 1-transitions are not possible for them. Indeed, the left branch 
of the parabola-like structure for the a-electrons (related to transitions where the £ quantum number decreases by 
one) shows gaps for the f = |m| orbitals whereas for the tt electrons it is the right branch where the corresponding 
lines are missing. The spectrum of the total dipole is shown above the contour plot. It is a broad structure fragmented 
in lines originating from single particle transitions. Around 20 eV there is a minimum in the total spectrum although 
the KSORD spectra show lines there. This is due to destructive interference. Instead, the dipole emission from KS 
orbitals with (initially) high £ quantum number at 17 and 24 eV interferes constructively, leading to two pronounced 
peaks in the total dipole spectrum. The small zero frequency component is caused by the few electron density which 
was freed by the delta-kick. 

The dipole spectrum in Fig.^ can be compared with the linear response results in |l^ (Fig- 3) and in |Q (Fig. 1). 
In Ref. |l4[ where a similar jellium model was studied within linear response theory we find similarities with our 
spectrum in the energy region 14-20 eV, in particular the pronounced 17 eV- line. Instead, in Fig. 3 of there is 
weak dipole emission below 10 eV (apart from a peak at 4eV) and also low dipole strength above 25 eV, apart from 
a broad but rather weak hump with a maximum around 34 eV. In our spectrum in Fig. ^ lines in the region 25-35 eV 
are much stronger because also transitions which are not directly accessible by single particle processes starting from 
the ground state configuration are allowed in the nonlinear treatment. Similar to our result, in Fig. la of Ref. a 
highly fragmented structure is visible in the particle-hole strength function for energies 4-15 eV. 

IV. EXCITATION AND IONIZATION 

We calculated the interaction of our jellium model with ten cycle sin^-shaped (with respect to the electric field) laser 
pulses in the frequency range from 6.8 up to 47.6 eV. The peak intensity was adjusted in such a way that intensity 
times pulse duration (energy per unit area) was held constant, i.e., in atomic units A^uji — 1.8375 x 10~^ where uji 
is the incident laser frequency and A is the vector potential amplitude. With such laser intensities the probability 
to remove the first electron remained below 10% for all frequencies. After the laser pulse was over we continued 
the propagation of all KS orbitals to allow for delayed ionization and free oscillations. The total simulation time 
corresponded to tend = 11-5 fs so that a real Cgo which has been ionized at most to C^q has no time for fragmentation. 
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The result is presented in Fig. |[ The contour plot shows the logarithmically scaled dipole strength vs. incident 
frequency uji and emitted frequency coe- Left to the contour plot the number of removed electrons 

iV,e,„ = 250 - / dVn(r,te„d) 

JGrid 

is plotted as a function of the incident frequency uji . In all cases delayed ionization was negligible compared to the 
electron density which was freed owing to the laser pulse. 

The maximum dipole strength is along the diagonal LUi = uje where the excitation was resonant. However, at 
LOi — uOf. = 20.5 eV there is a relatively weak dipole response but the ionization yield has a local maximum (dashed 
line). Instead, the local minima in ionization at 17 and 24 eV coincide with emission in two strong lines near the same 



values for We- Those lines were already observed in Section [II. For uji = 30 eV one observes both, strong ionization 
and a strong dipole response. For even higher uji ionization drops and the jellium system does not provide modes 
> 40 eV. 

The mutually exclusive behavior of ionization and dipole emission at 17, 20.5, and 24 eV might be surprising. One 
could expect that resonantly driven transitions yield big excursions of the electron density, thus leading to strong 
ionization also. This was observed for the case of sodium clusters interacting with strong laser pulses of 100 fs duration 
psf , and also the experimental results in [T^ ] were analyzed by unfolding excitation and (subsequent) autoionization. 
We believe that our different findings arc due to the short duration of the exciting laser pulse. Although at 17 
and 24 eV the short laser pulse excites resonantly Trcr-transitions with the outermost electron involved it does not 
efficiently couple those excitations to the continuum, thus keeping the ionization probability relatively low. Instead, 
at 20.5 eV the outer electron is not resonantly coupled to another bound state but rather absorbs a single photon to 
make a bound-free transition. The transitions that are resonant with the incident frequency of 20.5 eV involve low-£ 



a- and 7r-orbitals whose emitted radiation interfere destructively, as discussed already in Section [II . At 30 eV the 
outer electron can be freed by absorbing a single photon while stronger bound electrons are resonantly excited to 
perform bound-bound transitions of cr^-type which leads to efficient dipole emission also. Concluding this Section we 
can say that in very short exciting laser pulses the Cgo response seems dominated (i) by single particle transitions, 
although not exclusively those starting from the ground state, and (ii) by direct photoionization instead of ionization 
by resonantly driven single particle transitions or plasmons, as it is the case in longer pulses. 



V. IONIZATION MECHANISM AT A = 800 NANOMETER 



After studying the dipole response at relatively high energies we now turn to the interaction of Ceo with laser light 
of 800 nm wavelength. We simulated the interaction of our Cgo jellium model with a ten cycle 800 nm sin^-shaped 
laser pulses (corresponding to 26 fs pulse duration). In Fig. ^ the removed electron density A^rom after the pulse vs. 
the peak intensity of the pulse is presented. In the regime where A^rcm 1 this is equivalent to the single ionization 
probability. In the case of perturbative off-resonant multiphoton ionization of atoms it is well-known (see, e.g., [pTf) 
that the ionization probability is ~ /" where n is the number of photons needed to reach the continuum from the 
initial state. From Section ^ we know that in our jellium model the first ionization potential is Ip = 0.279, and 
thus we expect n = 5 photons necessary for ionization if Ceo behaves multiphoton-like. As is evident from Fig. ^ this 
is the case, in agreement with an experiment performed using a Ti Sapphire laser (A = 800 nm, 30 to 120 fs pulse 
duration) ||l^. There, from the Cqq yield-slope, also the second electron was found to behave multiphoton-like with 
n ^ 8. Those findings are in contrast with earlier experimental results in pT] (same wavelength and pulse durations) 
where excitation of the 20 eV Mie-resonance, corresponding to n = 13 photons, was concluded to be the dominant 
ionization mechanism. The n — 13 slope is depicted in Fig. ^ in the upper left corner and is much too steep to fit 
with our numerical result. This appears reasonable because the incident laser frequency 1.6 eV lies energetically far 
below 20 eV so that in short pulses an efficient excitation of the latter is unlikely. The experimental fragmentation 
onset and the Ceo saturation intensity as observed in [|o| is also indicated in Fig. ^. 

In Fig. ^ we present dipoles and the corresponding spectra for two particular laser intensities. For the higher 
intensity [plots (c) and (d)] the electron density continues oscillating quite strongly after the pulse but with little 
effect on ionization. In the spectra (b) and (d) the black curve was calculated from Fourier-transforming the dipole 
with respect to the entire time interval shown in the plots above. Therefore lines corresponding to the laser harmonics 
can be observed in those spectra. Fourier-transforming only over the time after the pulse leads to the gray curves. 
In the low intensity case (b) the dominant line is around 3eV, corresponding to the transition of the outermost 
TT-electron from ^ = 4 to ^ = 5. This transition is nearly resonant with two laser photons. In the high intensity case 
the electron density continues oscillating with the laser frequency even after the pulse. The next more energetic peaks 
are, again, the £ = 4^5 transition near 3eV followed by the £ = 9 ^ 10 transition of a cr-electron at 5.5 eV. The 
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next line corresponds to the ionization energy 7.5 eV. For both intensities excitations beyond 10 eV are several orders 
of magnitude weaker and can therefore not play any role in the ionization process of our Cgo model at that laser 
wavelength. 



VI. CONCLUSION 



In this paper we presented a nonlinear time-dependent density functional theory (TDDFT)-treatment of a Cgo 
jellium model in intense laser pulses. The Kohn-Sham (KS) orbitals were expanded in spherical harmonics and 
discretized in radial direction on a sufficiently big numerical grid on which they were propagated in time. By this 
method all bound states and the continuum can be properly represented, allowing for an accurate description of 
ionization and higher order transitions, not included in linear response theory (sometimes also called TDDFT). 

We characterized the KS orbital resolved dipole (KSORD) response. The KSORD spectra can be understood in 
terms of single particle transitions. The dipole radiation emitted by individual KS orbitals can interfere destructively. 
This is the case for tT7r-transitions between KS orbitals with low I (and £ ± 1). The dominant lines in the total dipole 
spectrum at 17 and 24 eV were found to originate from single particle tTTr-transitions with high-i* KS orbitals involved. 
We also found a strong response in the dipole spectra around 30 eV. This broad structure is fragmented in single 
particle transitions of (T(5-type. 

For short exciting pulses with frequencies between 7 and 48 eV we studied both ionization and the dipole spectra, 
calculated from the nonlinear dipole response after the pulse had passed by. We observed that when the incident 
frequency matched with the two resonant peaks of the dipole spectrum at 17 and 24 eV, respectively, ionization 
dropped. Instead, an incident frequency in between (w 20.5 eV) led to a maximum ion yield. This was attributed to 
the fact that in very short exciting laser pulses the Ceo response is dominated by direct photoionization instead of 
by ionization due to resonantly driven single particle transitions or plasmons. Hence, when the outermost electron 
is coupled resonantly to another bound state ionization decreases. A maximum in both ionization yield and dipole 
strength was found at frequencies around 31 eV where the outer tt electron is directly photoionized and stronger bound 
electrons are resonantly excited. For even higher incident photon energies > 37 eV ionization decreased because energy 
transfer to the Ceo became inefficient. 

For 800 nm fs laser pulses we observed direct multiphoton ionization rather than ionization via excitation of a "Mie- 
resonance" . This result agrees with recent experimental findings in . Although it disagrees with earlier results in 
PI we are quite confident that our model yields correct predictions in this respect because the multiphoton character 
for optical (or near optical) frequencies is well-known for metal clusters (see, e.g., and references therein) where 
collective electron behavior is more pronounced and resonances are less energetic. 

The TDDFT jellium model offers one of the very few feasible approaches for theoretical investigations concerning 
the interaction of intense laser light with complex systems where many electron-effects, bound-free transitions, and 
rescattering might be important. Preliminary results for "above threshold ionization" (ATI) of Ceo (as in the experi- 
ment 1^) were also obtained within this model [ p6| . Simulating ATI is always much more demanding than studying 
the dipole response and ionization only since one has to trace the free electron motion in the continuum as well as 
rescattering events with very high accuracy over several laser cycles of optical frequency. 
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TABLE I. The ground state single Kohn-Sham (KS) particle orbital energies as obtained on our numerical grid used for 
propagation (grid spacing Ar = 0.3). The a-electrons occupy values up to 9, the 7r-electrons up to 4. The 5-orbitals are not 
occupied in the ground state configuration. Orbital energies of unoccupied levels are written italic. The occupied orbitals are 
enumerated for the sake of easy reference. For each ^-shell £+1 KS orbitals to account for all the possible \m\ (numbered from 
lowest to highest) are needed. 



FIG. 1. 

Ground state poperties of the jellium Ceo-model. (a) The total potential V(r) + U[n{r)] + Vxc[n{r)] depends on 
the angular quantum number i (through the centrifugal barrier). For n = 1 (cr-electrons) orbitals from £ = up to 
£ — 9 are occupied in the ground state situation, for n = 2 (7r-electrons) orbitals from £ = to £ = 4 are occupied. 
The potential for the empty ^ = 10 orbitals is drawn dashed. The radial shape of the total ground state density n{r) 
is also plotted, (b) The single Kohn-Sham particle energy levels corresponding to the potentials in (a), cr-states are 

drawn solid while 7r-states are plotted dashed. The degeneracy is 2(2£+ 1). (c) The orbital densities n'f'^K The sum 
of those is ?i(r)/47r. 



FIG. 2. 

Kohn-Sham (KS) orbital resolved dipoles (KSORD, left) and the corresponding spectrum (right) after a delta-kick 
with an electric field E{t) = AS{t), A = 0.01. The Kohn-Sham orbitals are enumerated according Table |[ The total 
dipole d{t) and the total dipole strength (Dip.str) are given also. See text for a detailed discussion. 



FIG. 3. 

Right-hand-side: contour plot of the Fourier-transformed dipole d{t) after ten cycle pulses of frequency uji. The 
dipole strength vs. incident frequency uji and emitted frequency is logarithmically scaled (cf. color bar ranging 
from lO"'^" to 10~^, in arbitrary units). Left-hand-side: number of removed electrons. The horizontal solid (dashed) 
lines (line) indicate frequencies tOi where ionization was relatively low (high) and excitation was efficient (inefficient). 



FIG. 4. 

Removal of the first electron vs. the peak intensity of a ten cycle 800 nm sin^-pulse (solid curve). A multiphoton-like 
behavior ~ /" with n — 5 is evident. Instead, n = 13 photons would be necessary to excite a 20 eV Mie-resonance. 
Fragmentation threshold and C^g saturation intensity from [|lO| are also indicated (vertical lines) . 



FIG. 5. 

Dipole d{t) [plots (a) and (c)] and its Fourier transformation [plots (b) and (d)] for a ten cycle 800 nm sin^-pulse 
of peak intensity / = 4.56 x lO^^-yv/cm^ [plots (a) and (b)] and / = 2.85 x lO^^-v^/cm^ [plots (c) and (d)]. The 
black spectra are calculated from the entire time interval, the gray curves (multiplied by 10^, for better legibility) are 
spectra calculated from the time after the laser pulse. The former exhibit laser harmonics while in the latter single 
particle hole-excitations at low energies are dominant. Excitations > 15 eV are many orders of magnitude too weak 
for having a strong effect on ionization. 
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